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We present an experimental investigation of the non stationary frictional properties 
of multicontact interfaces between rough elastomers and rough hard glass at low 
velocities (< 200 /im s _1 ). These systems, for which the deformation contribution to 
friction is negligible, are shown to exhibit a phenomenology which is similar to what 
is observed for non elastomeric materials in the same multi-contact configuration, 
and which are quantitatively described by the state- and rate-dependent friction 
laws. This permits to identify clearly the two contributions to adhesive friction 
which are mixed in steady sliding : the interfacial shear stress which appears as 
thermally activated formation and breaking of molecular bonds, and the real area 
of contact which evolves through viscoelastic creep of the load bearing asperities. 
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1. Introduction 

Studies of friction of elastomeric materials over hard surfaces have been, up to now, 
limited to stationary sliding, although practical situations are often non stationary. 
Extensive experimental studies of the sliding friction of rubber show that the steady 
state dynamic friction force F ss varies with temperature T and velocity V according 
to a single 'master curve' (Grosch, 1963): 

F SS (V,T) = F ss {a T V) (1.1) 

where ax is a decreasing function of temperature, characteristic of the bulk vis- 
coelastic properties of a given elastomer (see eg. Ferry, 1970). Furthermore (Grosch, 
1963) for materials with widely different viscoelastic properties, this curve shows 
a single peak for smooth tracks (plane glass) with its maximum occurring at a 
velocity V s such that V s / fc (where fa" is the frequency at which the bulk loss 
modulus is maximum) defines a length d s of molecular size. For rough tracks, a 
second peak appears, with its maximum at a velocity Vh such that Vh/ fs (where fg 
is the frequency at which the bulk loss tangent is maximum) , defines a length dh of 
the order of the track roughness. Two different mechanisms were associated with 
these two length scales, one attributed to adhesion loss involving interfacial molecu- 
lar bonds, the other one to viscoelastic dissipation within the bulk of the asperities. 
For the adhesion component, the mechanism by which a frequency characteristic 
of the bulk might be involved is not clear, and several theories where proposed to 
explain the bell shaped master curve observed. 

f To whom correspondence should be addressed 
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Following the seminal adhesion theory of metal friction of Bowden and Tabor 
(1950), Ludema & Tabor (1966) described the adhesive friction force between a 
spherical hard indenter and a flat viscoelastic surface as a combination of two 
physical terms : 



where S r is the real area of contact and a s is a mean interfacial shear stress required 
to break junctions. In their single Hertz contact geometry, the area of contact varies 
E(to) 2 / 3 where E is the time dependent Young modulus evaluated at 
the time t necessary for the contact to travel its own size at velocity V, namely 
t = a/V where a is the contact size. As the viscoelastic modulus obeys a time- 
temperature equivalence with a shift factor ax and increases with reduced time 
t/a,T, S r increases with to far and decreases, in stationary sliding, with V. The 
shear stress a s is expected to increase with the reduced shear rate are as observed 
in bulk tensile fracture experiments on the same kind of elastomers (Smith 1958). 
Assuming that the shear is accommodated in a thin interfacial layer of thickness h, 
e = V/h, and a s increases with cltV . Finally, the friction force reads 



which produces a bell shaped curve with a position of the maximum in good agree- 
ment with experimental results, provided that h is chosen to be of order 10 nm 
(Ludema & Tabor, 1966). 

Savkoor (1987) used a similar approach to describe friction of rough surfaces, 
as the result of the statistical averaging, over the contact population, of a single 
micrometric contact response. He analyzed the two effects, with a more detailed de- 
scription of the evolution of the contact area, and introduced a fracture mechanics 
description of the shear strength, based on a Griffith-like energetic criterion for fail- 
ure of the contact. This theory, though leading to the correct qualitative behavior, 
depends highly on details such as geometric characteristics of the surfaces. 

Another approach, proposed by Schallamach (1953, 1963) focuses on the inter- 
facial term, which is described as thermally activated formation and breaking of 
bonds between rubber molecules and the hard substrate. This yields a mean force 
necessary to break bonds that increases with velocity, but a number of bonds that 
decreases with velocity, leading to a global friction force that qualitatively has the 
observed bell shape. 

In all these approaches, adhesive friction is viewed as a competition between 
two effects : the material strength (force to break a bond, shear strength) increasing 
with velocity, and the amount of adhesive (number of bonds, area of contact) that 
decreases with velocity. In steady state, these effects are mixed and thus difficult 
to separate. 

Recently, studies of non stationary friction of non elastomeric materials have 
been successful in separating these two contributions to the friction force. They 
concern friction at low velocities (< lOO^m s _1 ) between macroscopically flat sur- 
faces with a micrometric roughness, forming a multicontact interface (MCI) with 
a large but dilute population of contacts. The phenomenology of these systems, 
which has proved to be independent to a large extent of the nature of the materials 
involved (granite, paper, polymer glasses), is quantitatively described by the state- 
and rate- dependent friction (SRF) laws first introduced by Rice & Ruina (1983). 



F = u s E r 



(1.2) 



F ss (a T V) = <j(a T V)Z(to/a T = a/{a T V)) 



(1.3) 
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The detailed analysis of polymer glasses has shown that these laws can be inter- 
preted through the Tabor-like decomposition (1.2) involving the real area of contact 
S r and the interfacial shear strength a s . For multicontact interfaces, this decompo- 
sition is particularly well suited to explain the Amontons- Coulomb proportionality 
observed between the friction force and the normal load through a linear dependence 
of the real area of contact on normal load, coming from the statistical distribution 
of asperity heights (Greenwood & Williamson, 1966). The small variations of the 
friction coefficient which are responsible for the dynamical stability of the system 
depicted on figure 1 were extensively studied. They could be attributed to the evo- 
lution of the real area of contact with a "state" variable — a time-dependent age <f> 
— and to the dependence of the interfacial shear strength a s on the instantaneous 
shear ratef, which is itself proportional to the instantaneous velocity x of the slider. 
That is, the friction force can be written : 

F(x, </>) = WfM(x, 0) - tr 8 (i)E r (0) (1.4) 

The variations of the contact area with age could be obtained by measuring the 
increase in static friction with the time of contact t c prior to sliding. For polymer 
glasses, this is associated to the plastic creep of the load-bearing asperities, as 
observed by direct visualization of the real area of contact (Dieterich & Kilgorc, 
1994) and by quantitative comparison with the bulk plastic creep properties for 
various temperatures below the glass transition (Berthoud at at, 1999). In steady 
sliding, this ageing occurs during a time Dq/V where Do is a characteristic length 
found to be of the same order of magnitude as the mean contact size. The friction 
force then has a form similar to Ludema & Tabor's : 



F SS (V) = o- s (V)X r (D Q /V) (1.5) 

which, at low velocity (< 100 /xms -1 ), was found to be a decreasing function 
of velocity, leading to the existence of a bifurcation between steady sliding and 
stick-slip oscillations. This bifurcation could be analyzed quantitatively with a time 
evolution law for the state variable <f> in non steady conditions that accounts for 
the past history of the system over the length D , and that interpolates between 
the static and steady state regimes : 

« t) _£^(_£»zi<Q) « (1 . 6) 

which takes the differential form : 

The dependences of real area of contact on age and of shear strength on velocity 
were found to be logarithmic : 

<j s (x) = a° s (l + a\n(x/Vo)) (1.8) 
E r (0) = £?(l + mln(^M>)) (1-9) 

f provided that no mechanisms of strength evolution with time, such as polymer chains inter- 
diffusion, is at work (Berthoud et al, 1999) 
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with values of a and m small compared to unity. The friction coefficient thus re- 
duces, to first order in a and m, to the usual SRF expression : 



When the steady state friction decreases with velocity (i.e. when B — A > 0), this 
friction laws predicts (see Appendix A) a bifurcation from steady sliding to stick- 
slip oscillations at low system stiffness and/or high normal loads for the system 
depicted on figure 1, as observed experimentally. The detailed study of this bifur- 
cation, coupled with steady state and static measurements, gives access to all the 
parameters involved in the SRF laws. This kind of analysis, performed as a function 
of temperature for polymer glasses, showed that the shear response of the contacts, 
which is expected to occur at the junction between asperities within a molecular 
thick layer of weaker material (confined polymer chain ends) is an activated process 
with an activation volume of order (1 nm) 3 , indicating that this nanometer thick 
layer between the contacts slides by activated depinning of uncorrelated nanomet- 
ric blocks, giving a picture very similar to Schallamach's description of interfacial 
sliding of rubbers (Schallamach, 1953, 1963). 

This paper deals with the extension of this description of non stationary friction 
to elastomers. We show that the low velocity friction properties of a multicontact 
interface between elastomers and rough hard glass, studied as a function of tem- 
perature, are indeed quantitatively described by the SRF laws. The viscoelastic 
properties of the material appear through a temperature dependence of the SRF 
parameter A and B compatible with a WLF transform. This permits to separate 
clearly the shear strength and the area of contact contributions to the friction force, 
and to bridge between Tabor and Schallamach's descriptions of elastomer friction. 
The paper is organized as follows : in the next section, we describe the experimental 
setup. Section 3 presents the results of the static, steady state sliding and unsteady 
friction measurements, followed by a discussion of these results within the SRF 
framework. 



We studied three elastomers made of synthetic poly-isoprene, reticulated with 
dicumyl peroxide at 443 K. The glass transition temperature T g , determined by 
differential scanning calorimetry, was changed by adding plasticiscrs (liquid paraffin 
or resin). The samples were moulded in 2.8 mm thick plates. The composition in 
mass and the glass transition temperatures of the samples are shown on table 1. 

The dynamic mechanical properties (storage G and loss moduli G") of the three 
elastomer samples were obtained in the frequency range 0.1 — 750 Hz for various 
temperatures. This permitted to perform a WLF transform (shown on figure 2(a) 
for sample #1), leading to the shift factors ax that are presented, on figure 2 with 
the glass transition as a reference temperature. Most of the results obtained are, 
unless stated, qualitatively independent of the sample and we will therefore present 
the quantitative results obtained on sample #1. 




(1.10) 



2. Experimental 



(a) Samples 
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Table 1. Sample composition and glass transition temperature 





Sample #1 


Sample #2 


Sample #3 


poly-isoprene 


100 


100 


100 


paraffin 


20 






resin 




27 




peroxide 


1.2 


1.2 


1.2 


glass transition (K) 


303 


296 


274 



The slider is composed of an aluminum plate (lOx lOx 1 cm 3 ) on which was 
glued the elastomer sample (7 x 7 x 0.28 cm 3 ). The surface of the slider was hand 
lapped at low temperature (below T g ) with a SiC abrasive powder of nominal grain 
size 23 ^m. This procedure has proved to give a reproducible rms roughness of 
1.3 fim. 

The track (15 x 15 x 1 cm 3 ) is made of rough glass prepared in the same way. 
We use a glass track in order to avoid any dependence of the shear strength a s with 
the contact age 4> that might be caused by inter-diffusion of polymer chains. 

(b) Friction apparatus 

The details of the experimental setup sketched on figure 1 are described in 
detail elsewhere (Berthoud et al. 1999). The slider is driven along the track by 
a translation stage, composed of a stepping motor pushing a cantilever spring of 
stiffness k. The spring pushes the slider through a cross-cylinder junction in such 
way that the contact is in the same plane as the track-slider interface. The velocity 
range of the driving stage is V = 0.1 to 200 (im s^ 1 ensuring that frictional heating 
is negligible. The spring stiffness k was chosen to be the most compliant part of 
the setup in order to be a control parameter of the stability of the system. Using 
the relaxed storage modulus of the most compliant sample, we can estimate a lower 
limit for the shear stiffness of the samples in the temperature range studied of 10 6 
N m _1 . We used cantilever springs of maximal stiffness 2.4 x 10 5 N m -1 . The force 
applied to the slider is deduced from the spring deflection measured by an eddy 
current inductive transducer with a sub-micron resolution within the frequency 
range. 

The temperature control of the samples was achieved by gluing resistive stripes 
(total of 140 W) under the track and over the slider, each with a 100S1 platinum 
probe and a controller for thermal regulation. The elastomer temperature was also 
measured with a separate probe located inside it. Both track and sample where 
enclosed in a styrofoam box. The temperature was thus controlled between room 
temperature and 370 K with 0.5 K precision. It was found that humidity has an 
influence on the results only when above 60% and the results presented here were 
obtained under dryer conditions. 

(c) Friction force measurement 

This setup allows the measurement of both static and dynamic friction forces. 
The static force F s is defined from the peak value of the tangential force when 
loading rapidly the slider from rest. Both static and dynamic friction forces show 
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similar fluctuations that depend on the slider's position along the track and which 
correspond to long wavelength inhomogeneities of the interface. They can thus be 
removed by measuring the friction forces AF relative to dynamic friction in steady 
state at a reference velocity V If} { measured at the same position. This enabled us 
to measure the variations of friction forces with an error lower than the track 
fluctuations. The absolute value of the force is obtained by adding to AF the mean 
value of the friction force at the reference velocity V lc f over the whole track. 

Figure 3 shows the static and dynamic friction forces measured as a function of 
the normal load W. 

The proportionality defines the static and dynamic friction coefficients as the 
slopes fi s = F s /W and [id = F4/W. This proportionality of the friction force to the 
normal load is usually attributed to a linear relation between the real area of contact 
and the normal load, as expected from Greenwood's description of multicontact 
interfaces (Greenwood & Williamson, 1966). 

3. Results 

(a) Static friction 

In order to study the dependence of the static friction force on contact time t c 
prior to sliding, we follow a procedure proposed by Bcrthoud et al. (1999) : the slider 
is first put in steady sliding at a velocity V YC f during 10 seconds in order to renew 
the contact population. It is then stopped, remaining under tangential load for a 
resting time t c . After that time, the slider is shear- loaded at velocity Vj oa d and the 
force is recorded. The static friction is then defined from the peak value of the force, 
at which point the instantaneous velocity of the slider equals the loading velocity. 
It has been shown (Berthoud et al. 1999) that the loading rate has an effect on the 
value of the static friction coefficient. We thus measure it at the constant loading 
rate Vi oa d = 100/im s . For convenience, the reference velocity V le { is chosen equal 
to the loading velocity. Figure 4 shows the static friction coefficient function 
of contact time t c for different temperatures. We find a quasi-logarithmic increase 
of the static friction with rest time, with a mean slope depending strongly on 
temperature. In order to analyze this dependence within the SRF framework, we 
define a local logarithmic slope (5 S = dfi s /dlnt c over the first two decades of time. 
Figure 5 shows the dependence of (3 S on T, showing a clear maximum around 
60° C. In order to analyse the influence of temperature on this ageing effect, we 
tried to check whether it obeys a time-temperature equivalence. By applying a 
temperature-dependent scaling 1/ay in time to each set of measurements, we could 
obtain a single 'master curve' (figure 6(a)) for the static friction coefficient : 

H s (t c ,T) = n s {t c /a s T ) (3.1) 

Furthermore, figure 6(b) shows that the shift factors af. agree quantitatively with 
the viscoelastic bulk ones. Static ageing is thus the result of an increase of the real 
area of contact through viscoelastic creep of the load bearing asperities. 

(b) Steady state dynamic friction 

The steady state dynamic friction coefficient /xjf(V) was measured for various 
temperatures, as shown on figure 7. The friction always decreases with velocity 
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within the temperature range studied, but with slopes that depend on temperature. 
Again, we define a local logarithmic slope (3d — —d/j, d s /dlnV between 0.1 and 10 
/im s _1 , which is shown as a function of temperature on figure 8. It exhibits a clear 
maximum at about 330 K and the trend at higher temperatures indicates that (3d 
might become negative. As for static friction, we are able to define temperature- 
dependent scaling factors in time, with the help of which measurements at 
various temperatures fall on a single 'master curve' shown on figure 9(a) : 



Figure 9(b) shows that the scaling factors compare quantitatively with bulk ones, 
in agreement with previous sliding friction measurements on various elastomers 
(Grosch, Ludema & Tabor), although we don't see the peak that was observed in 
these experiments. However, the saturation found at low reduced velocity indicates 
that such a peak might be present in our system at higher temperatures. 



Because of the velocity weakening behavior of the steady friction force, station- 
ary sliding could only be obtained under low loads and with stiff external springs. 
The system otherwise slides through a succession of sticking and slipping phases. 
The bifurcation between these two dynamical regimes is continuous : the ampli- 
tude of the stick-slip oscillations continuously grows from zero when increasing the 
normal load as shown on figure 10 (a). 

Figure 10 (b) shows the dynamical phase diagram of the system in the (k/W, V) 
space, at room temperature. The critical reduced stiffness slightly depends on the 
velocity V as observed on other materials in the same multicontact configuration. 

This kind of bifurcation is accounted for by the linear stability analysis of the 
SRF laws (see Appendix A). It is predicted to occur for 



The critical reduced stiffness (k/W) c was measured at a driving velocity V = 
2/ims -1 for various temperatures, together with the pulsation f2 c of the stick-slip 
oscillations at the bifurcation. Using the previous measurements of the velocity 
weakening slope (3d at corresponding temperatures, we can then deduce the value 
of the memory length D = (3d/ (k/W) c . 

Figure 11 shows the values obtained for different temperatures for the three 
elastomers used. Though the dispersion is large due to cumulative errors on (3d and 
{k/W) c , the memory length appears independent of temperature and sample type, 
having an average value of D ~ 1 fim, comparable with the surface roughness, 
in accordance with the geometrical interpretation of this length in terms of a slip 
distance necessary to renew the contact population. 

The pulsation O c gives access to the dependence of the friction coefficient on 
the instantaneous velocity A = d[id/d\nx as (see Appendix A) 



(3.2) 



(c) Non stationary friction, stick-slip bifurcation 




(3.3) 




(3.4) 
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Its temperature dependence is shown on figure 12. The indirectness of the mea- 
surement leads to cumulative errors that become large, but there is nevertheless a 
clear increase with temperature. 



4. Discussion 

We have so far shown that the velocity weakening regime in steady sliding lead- 
ing to a continuous bifurcation between steady sliding and stick-slip oscillations is 
compatible with the phcnomcnological state- and rate-dependent friction laws. This 
permitted the determination of the various parameters involved. A further step can 
be made by trying to link the dynamic response to the static threshold. Provided 
that the loading phase is short compared to the static age t c of the multicontact 
interface, static friction is a transient measurement at the loading velocity Vi oa d 
with an age t c , and following the SRF description, the static friction coefficient can 
be written as 

Hs{tc) =^o + A\n{Vx oaA /V )+B\n{t c /t ) (4.1) 

This implies a logarithmic dependence of the static friction on the loading veloc- 
ity, which was indeed observed and found quantitatively compatible with this law 
for polymer glasses (Berthoud et al. 1999). Furthermore, the logarithmic ageing 
slope (3 S should compare with B. As we have measured independently the direct 
effect parameter A, and the steady state friction slope (3d = B — A, the comparison 
can be done as shown on figure 13 where (3d + A and (3 S are plotted for different 
temperatures. The agreement indicates that the state-and rate-dependent friction 
model coherently describes the friction of rough elastomers on rough glass. Further- 
more, the memory length D was found to be independent of either temperature or 
material and to be on the order of the mean contact size. This shows that the inter- 
pretation of the age parameter in terms of real area of contact which was shown to 
hold for polymer glasses also holds here, though the physical origin of the ageing is 
clearly viscoelastic in our case, as was suggested by Ludema & Tabor (1966). This is 
confirmed by the fact that the dependence of static friction on both time and tem- 
perature follows a time-temperature equivalence compatible with bulk viscoelastic 
measurements . 

The comparison of our results with other experimental studies of steady sliding 
friction on elastomers is limited because, with our system, we only observed velocity 
weakening behavior. However, we have seen that the saturation of fi s £ at low reduced 
velocities (figure 9(a)) may be the trace of the maximum observed in previous 
investigations. We assumed that the viscoelastic dissipation due to the deformation 
of asperities while sliding was negligible compared to adhesive dissipation. This 
can be checked by an evaluation of the deformation contribution to the friction 
coefficient as a function of viscoelastic properties of the material and the geometric 
characteristics of the surfaces forming the multicontact interface. At sliding velocity 
V, the contacts are stressed during a time of the order of V/a where a is the mean 
contact size. Such a contact will be subjected to a deformation e ~ a/R, where R is 
the asperity radius of curvature, within a volume of order a 3 . The energy dissipated 
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by viscoelastic losses in the contact will then be of order 

E dei ~ G"(V/a) x e 2 x a 3 ~ G"(V/a)j^ (4.2) 

This dissipation will result in a contribution F de f to the friction force such that 
afdcf = ^def, or to a contact shear resistance 

a dcf ~ G"-^ (4.3) 

From Greenwood's (Greenwood & Williamson, 1966) description of elastic multi- 
contact interfaces, we can estimate the mean contact size as a ~ \J Ra and the 
mean contact pressure as 

P = W/X r = G'(^) 1/2 (4.4) 

where a is the characteristic length of the asperity height distribution. The de- 
formation loss will thus give a contribution to the friction coefficient of the order 
of 

li dei {a T V) ~ ~ (^) 7 (tan<5) aT/=y/v ^ (4.5) 

This is in agreement with the experiments of Grosch (1966) on dry silicon carbide 
paper for which a and R are expected to be of the same order of magnitude. Note 
that a similar result was obtained by Persson (1998) but in the limited case of very 
rough surfaces for which a and R arc of the same order of magnitude. For our 
surfaces, these quantities, measured from surface profiles, are of order a ~ 1 /im 
and R ~ 100 fim. Furthermore, the loss tangent of sample #1 has a maximum of 
0.9 about a reduced frequency arfs — 10~ 3 Hz. The deformation component of 
friction is thus expected to have a maximum of order 0.09 at a reduced velocity 
arV ~ axfsVo-R ~ 10~ 2 ^m s _1 . From figure 9, we see that there is no peak about 
this velocity and that the measured friction is about five times higher. This is also 
the case for the other samples and the observed friction comes clearly from adhe- 
sive mechanisms in the contact junctions. From our interpretation of this adhesive 
component of friction, the change from a regime of velocity weakening to velocity 
strengthening at high temperatures or low velocities would arise from an inefficient 
ageing (state) effect compared to the rheological shear response (rate) of the inter- 
face. This gives further support to the adhesion view of Ludema & Tabor. We thus 
propose an extension of their description of elastomer friction accounting for non 
steady dynamics as 

F(x,<P,T) =<7 s (a T ±)£r( — ) ( 4 -6) 

(it 

with the time evolution of the state variable <j> given by equation (1.7), which 
reduces to equation (1.3) in steady state, the contact size being averaged over the 
contact population. We don't expect the interfacial shear strength to compare with 
the bulk shear strength of the elastomer as no wear effect could be observed. A 
mechanism similar to Schallamach's (1953) description of interfacial friction, in 
terms of thermally activated and uncorrelated depinning of interfacial molecular 
bonds, as was proposed for glassy polymers, is more probable. 
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5. Conclusion 

We have shown that the non steady adhesive friction of elastomers over hard glass, 
with a well defined surface geometry forming a multicontact interface is quantita- 
tively described by the phenomenological state- and rate-dependent friction laws. 
This permits the separation between the two contributions to the friction force : 
the interfacial shear strength and the area of contact, which arc mixed in steady 
sliding. This extends Ludcma & Tabor's (1966) description of elastomers friction in 
terms of an interfacial shear strength similar to that invoked by Schallamach (1953, 
1963), and a real area of contact that was shown to evolve through viscoelastic 
creep within the bulk of the contacting asperities leading to memory effects over a 
length characteristic of the mean contact size. 

We are grateful to T. Baumberger and Ch. Caroli for fruitful discussions and critical 
reading of the manuscript. The Groupe de Physique des Solides is "associe au Centre de 
la Recherche Scientifique et aux Universites Paris 6 et Paris 7" . 



Appendix A. 

This appendix presents the linear stability analysis of steady sliding at velocity V 
of the dynamical system shown on figure 1 with a friction force between the slider 
and the track described by a friction coefficient depending on the instantaneous 
velocity x of the slider, and on the state variable 0, the time evolution of which is 
given by equation (1.7). At the low velocities considered, inertia can be neglected 
and the dynamical equations of the system are : 

-W/j, d (x,(j>) + k(xo + Vt-x) = (Al) 

= 1-^ (A2) 



where xo is the spring elongation in steady sliding at the driving velocity, for which 

(A3) 
(A 4) 



x = V (A3) 

= ^ 

v V 

^(V) = ii d (V,D /V) (A 5) 



We are looking for the dynamical evolution of small perturbations about the steady 
state solution 

x = V + Sx (A 6) 

cj> = y + H (A 7) 

with 5x <C V and 8<p <C D /V Substituting into the dynamical equations, and 
expanding to first order about the steady state solution, we find 

i)T V 

f*=-- (A9) 
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where partial derivatives are evaluated at the steady state value (x = V, <fr = D /V). 
Looking for solutions of the form 

Sx = Sx e st (A 10) 

5<fr = 6<fr e st (All) 

leads to a linear system of two equations in 5xq and 5 (fro that has non trivial solutions 
if 



2 d ^d . 

s + s 
ox 



where we have used the relation 

<ws_ = 

dV 



k V d^ d s 
W + Do~~dV~ 



+ 



kV 
DoW 



dju _ Do d^ d 

dx V 2 dcfr 



x=V 
<P=D /V 



Introducing 



B — A 
A 



d\nV 
djJLd 



dlnx 



x=V 
4>=D /V 



(A 12) 
(A 13) 

(A 14) 
(A 15) 



which might here depend on V, whereas in the usual SRF model are assumed to 
be constant, equation (A 12) reads 



2^ 

s v + s 



k_ 

W 



B-A 



Do 



kV 

+ DoW=° 



(A 16) 



The stability of the solutions (A 10) are given by the sign of the real part of the 
roots of this equation. A positive real part leads to an amplification of the pertur- 
bation an to instability of the steady sliding solution, whereas a negative real part 
implies the attenuation of the perturbation and steady sliding is stable. If A < 0, 
equation (A 16) has two real roots of opposite sign : 



s± = - 



V 
2A 



k_ 

W 



B-A 



Do 



± 



2T 



k_ 

W 



B-A 



-4. 



Ak 
'D W 



(A 17) 



the positive one leading to an unstable solution (A 10), and steady sliding is always 
unstable. If A > 0, the two roots are complex conjugate : 



s± = - 



V_ 

2A 



k_ 

W 



B-A 



Do 



± 



iV 



Ak 



2A V DoW 



k_ 

W 



B-A 



(A 18) 



a bifurcation happens when the real part of s± changes sign, steady sliding is stable 
if 



W 



> 



1 

D~o 



(B-A) 



(A 19) 
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and unstable otherwise. Furthermore, the unstable solution is oscillating at a pul- 
sation given by the imaginary part of the roots s±, which reads, at the bifurcation : 
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Figure 1. Sketch of the experimental setup : multicontact interface between nominally 
flat surfaces with a micrometric roughness. Slider = elastomer, track = glass. The slider, 
under normal load W is pulled by a spring of stiffness k, whose end is driven at constant 
velocity V. 

Figure 2. (a) Storage G' and loss moduli G" as a function of reduced frequency Ot/ for 
sample #1 (T g = 303 K). (b) Shift factors ar with the glass transition as a reference 
temperature for the three samples (□ sample #1; A sample #2; o sample #3). 

Figure 3. Static friction force (•) after a resting time of t c = 120 s, and steady state 
dynamic friction force (o) at velocity V = 100 /jm s~\ as a function of the normal load W 
for sample #1 at room temperature. The constant slope defines the static and dynamic 
friction coefficients fj, s and /xj. 

Figure 4. Static friction coefficient of sample #1 as a function of the resting time t c for 
three temperatures (o 293 K; x 333 K; ■ 368 K). At not too large times, the static ageing 
effect is clearly logarithmic, and the slope depends strongly on temperature. 

Figure 5. Logarithmic slope [3 S of the static ageing in the range 10 - 10 3 s, as a function 
of temperature relative to the glass transition temperature (sample #1, T g = 303 K). 

Figure 6. (a) Static friction coefficient of sample #1 as a function of the reduced contact 
time tc/oT, with scaling factors a a T computed to obtain a single curve (o 293 K; □ 313 K; 
o 323 K; x 333 K; + 338 K; A 343 K; • 353 K; ■ 368 K). (b) Comparison between these 
shift factors (•) and those obtained from bulk viscoelastic measurements (x). 

Figure 7. Dynamic friction coefficient fj, s d s as a function of the steady velocity V for three 
different temperatures : o 293 K; o 333 K; + 369 K. 

Figure 8. Local logarithmic slope (3d of the velocity dependence of the friction coefficient 

as a function of temperature. 

Figure 9. (a) Steady state dynamic friction coefficient as a function of the reduced ve- 
locity a^V (reference temperature T g ). The figure shows the measurements at different 
temperatures o 293 K; □ 308 K; o 333 K; x 353 K; + 368 K; (b) Shift factors from (a) • 
compared with bulk ones x (figure 2). 

Figure 10. (a) Friction force as a function of time for a driving velocity of V = 2 fim s _1 
for different values of the normal load. The steady state sliding is observed at low loads, 
whereas stick-slip oscillations occur at higher loads, with an amplitude that continuously 
increase from zero above the critical load, (b) stability diagram of the system as a function 
of driving velocity V and stiffness over normal load k/W. 

Figure 11. Memory length Do measured from the steady-sliding — stick-slip bifurcation 
at different temperatures for the three elastomer samples (• sample #1; □ sample #2; o 
sample #3). 

Figure 12. Direct effect parameter A measured from the bifurcation at V = 2/im.s" 1 as a 

function of temperature. 
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Figure 13. Comparison between independent measurements of the ageing component of 
the friction : from static friction f3 s (x) and dynamic friction (3d + A (o). 
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